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Abstract 

We perform evolutionary calculations of binary stars to find progenitors of systems with pa- 
rameters similar to the recurrent nova U Sco. We show that a U Sco-type system may be formed 
starting with an initial binary system which has a low-mass carbon-oxygen white dwarf as an ac- 
cretor. Since the evolutionary stage of the secondary is not well known, we calculate sequences 
with hydrogen-rich and helium-rich secondaries. The evolution of the binary may be divided into 
several observable stages as: classical nova, supersoft X-ray source with stable hydrogen shell 
burning, or strong wind phase. It culminates in the formation of a massive white dwarf near the 
Chandrasekhar mass limit. We follow the chemical evolution of the secondary as well as of the 
matter lost from the system, and we show that observed '^C/'^C and N/C ratios may give some 
information about the nature of the binary. 

1 Introduction 

Recurrent novae are a small class of objects which bear many similarities to other cataclysmic variable 
systems. They experience recurrent outbursts at intervals of 10-80 yrs. 

Recurrent novae have been proposed as possible progenitors of type la supernovae (Starrfield, 
Sparks & Truran 1985). Although they experience recurrent outbursts it is thought that the white 
dwarf mass would grow to finally exceed Chandrasekhar mass limit. 

Webbink et al. (1987) discussed the nature of the recurrent novae, and they concluded that accord- 
ing to outburst mechanisms there are two subclasses of these systems: (a) powered by thermonuclear 
runaway on the surface of the white dwarf (e.g. U Sco), and (b) powered by the transfer of a burst of 
matter from the red giant to the main-sequence companion. 

U Sco is one of the best observed recurrent novae. Historically, its outbursts were observed in 
1863, 1906, 1917, 1936, 1979, 1987 and in 1999. Schaefer (2004) discovered on Harvard College 
Observatory archival photographs unknown eruption of the recurrent nova U Sco in 1917. Schaefer 
notes that U Sco has a fairly constant recurrent cycle of 8-12 yrs, with about 25% of the outbursts 
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being missed due to proximity to the sun (including potential missed outbursts around 1926 and 1957). 
Determinations of the system visual luminosity at maximum and minimum indicate a range Amy ~ 9. 
Schaefer (1990) and Schaefer & Ringwald (1995) observed eclipses of U Sco in the quiescent phase, 
and determined the orbital period Porb=l -23056 d. 

For the 1979 outburst, ejecta abundances have been estimated from optical and UV studies by 
Williams et al. (1981) and Barlow et al. (1981). They derived extremely helium rich ejecta He/H~2 
(by number), while the CNO abundance was solar with an enhanced N/C ratio. From analysis of the 
1999 outburst, Anupama & Dewangan (2000) obtained an average helium abundance by number of 
He/H~0.4±0.06, while lijima (2002) 0.16±0.02. The estimated mass of the ejected shell for 1979 
and 1999 outbursts is ~ 10"'' M© (Williams et al. 1981; Anupama & Dewangan 2000). Hachisu 
et al. (2000) and Matsumoto, Kato & Hachisu (2003) estimated the envelope mass at the optical 
maximum as being ~ 3 x 10~^ M©. If this is the present case, the mass accretion rate should be 
smaller than ~ 2.5 x 10~^ M© yr~^ in the quiescence between 1987 and 1999. Spectroscopically, U 
Sco shows very high ejection velocities of (7.5-ll)x 10^ km s~^ (Williams et al. 1981; Munari et al. 
2000). Latest determinations of the spectral type of the secondary indicate a K2 subgiant (Anupama 
& Dewangan 2000; Kahabka et al. 1999). According to Kahabka et al. (1999) the distance to U Sco 
is about 14 kpc. 

Recently, the spectroscopic orbit of U Sco was determined by Thoroughgood et al. (2001). The 
radial velocity semiamplitudes for primary and secondary stars yield M^j = 1.55 ± 0.24 M© for the 
white dwarf and M2 = 0.88 ± 0. HM© for a secondary star. 

According to a model proposed by Kato (1996), supersoft X-ray emission should be observed 
about 10-60 days after the optical outburst. Indeed, BeppoSAX detected X-ray emission from U Sco 
at 0.2-20 keV just 19-20 days after the peak of its optical outburst in February 1999 (Kahabka et al. 
1999). The fact that U Sco was detected as a supersoft X-ray source (SSS) is consistent with steady 
hydrogen shell burning on the surface of its white dwarf component. 

In this paper we construct a grid of evolutionary sequences which may lead to formation of U 
Sco-like systems. In Section 2 we discuss the major phases of semidetached binary evolution. In 
Section 3 the evolutionary code is briefly described. Section 4 contains the results of the calculations. 
A general discussion and conclusions follow. 



2 The major modes of semidetached binary evolution 

We can identify three different modes of semidetached evolution of a close binary system with white 
dwarf as accretor: 

i) novae (nova outbursts); 

ii) stable hydrogen shell burning on a compact white dwarf; and 

iii) a strong, optically thick wind powered by steady hydrogen burning. 

While calculating evolutionary models of binary stars, we must take into account mass transfer 
and associated physical mechanisms which lead to mass and angular momentum loss. We can express 
the change in the total orbital angular momentum (7) of a binary system as 
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where the terms on the right hand side are due to: magnetic stellar wind braking (MSW); nova outburst 
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angular momentum loss, which describes the loss of angular momentum from the system due to non- 
conservative evolution (NOAML); optically thick wind angular momentum loss from the white dwarf 
during stable hydrogen shell burning (FWEND); and factional angular momentum loss in a strong 
optically thick wind (FAML). 

To account for magnetic stellar wind braking, we use the Skumanich (1972) 'law'. 



- 1 (5.0 X 10-2'') (271)^ G"^} 

where M^d and Mtot denote, respectively, the mass of the white dwarf primary and the total mass 
of the system; R^g is the radius of the secondary star; k2 is the radius of gyration of the secondary 
star, ^2 = 0-1 (Webbink 1976) and X characterizes the effectiveness of dynamo action in the stellar 
convective zone. Values of X = 0.73 (Skumanich 1972) or X = 1.78 (Smith 1979) can be found in the 
literature; however, we use a value of X = 1 according to Ergma & Sarna (2000) model calculations 
(the best fit to the observed orbital decay of the PSR B1744-24A). We assume that MSW losses are 
present in all three modes of evolution identified here. 



MSW 



2.1 Novae 

In novae, angular momentum loss accompanies mass loss due to the nova outbursts. 
2.1.1 Angular momentum loss associated with nova outbursts 

To take account of the angular momentum loss accompanying mass loss in nova outbursts, we use 
a formula based on that used to calculate angular momentum loss via a stellar wind (Paczyiiski & 
Ziolkowski 1967; Ziolkowski 1985 and De Greve 1993), 



= -^i-^2TfV^, Where M = f,M,g. (3) 

NOAML 'VHg'yhoi 
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Here, f\ is the ratio of the mass ejected by the white dwarf to that accreted by the white dwarf. 
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where AMacc and AMgj are the amounts of matter transferred from the secondary to the white dwarf 
and ejected in a nova outburst, during one nova cycle, respectively. The parameter /2 described the 
effectiveness of angular momentum loss during mass transfer (Sarna & De Greve 1994, 1996); Mgg 
denotes the mass of the subgiant secondary; M^g is the rate of mass loss from the secondary (—Msg 
is equivalent to the mass transfer rate). We employ the code developed by Marks & Sarna (1998) 
(hereafter MS98), which utilizes the results of the theoretical novae calculations of Prialnik & Kovetz 
(1995) and Kovetz & Prialnik (1997) (hereafter PK95 and KP97). Using this data, we are able to 
compute the value of f\, which we have previously had to assume as free parameter (Sarna, Marks 
& Smith 1996). By interpolation from the data of PK95 and KP97, at each time step, we use the 
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white dwarf mass and the mass transfer rate to determine the nova characteristics: /i, amplitude of 
the outburst (A), recurrence period (Pkc), and chemical composition of the ejected material (for more 
details see MS98). 

We take /2 = 1.0, which is typical for a stellar wind. Livio & Pringle (1998) proposed a model 
in which the accreted angular momentum is removed from the system during nova outbursts, as pre- 
viously suggested by Marks, Sarna & Prialnik (1997). Equation (3) is in quantitative agreement with 
these estimations. 

As discussed by MS98, angular momentum loss rates from nova outbursts are comparable with 
MSW losses. As an illustration, we can estimate the characteristic time scales [x = {j/J)~^] for both 
mechanisms. We assume a binary system in the same evolutionary stage as U Sco, with components 
of 0.94 and 1.3 Mq and orbital period 1.23 d. The secondary has a radius 2.22 Rq; fi — 0.92 and 
the mass transfer rate from the secondary to white dwarf is 3.3 x 10~^ M© yr~^ We obtain Xmsw 
1 .5 X 10^^ s and Tnoaml ~ 1 .7 x 10^^ s, i.e., angular momentum losses associated with nova outbursts 
is about an order of magnitude greater than those from the magnetic stellar wind. 

2.1.2 Accretion of material ejected during nova outbursts 

To calculate the re-accretion by the secondary of material ejected during nova outbursts, we assume 
that the mass of the material re-accreted (Mre-acc) is proportional to the mass of the material ejected 
by the white dwarf, such that 



The constant of proportionality is the ratio of the cross-section area of the secondary star to the area of 

a sphere at radius a from the white dwarf. We base this formula on the assumption that nova ejections 
are spherically symmetric, with ejection velocities large compared with the orbital velocity of the 
companion. 

During nova outbursts, matter which is re-accreted by the secondary changes significantly its 
chemical composition. 

2.2 Stable hydrogen shell burning 

We assume that during stable hydrogen shell burning, only angular momentum losses due to the 
magnetic stellar wind will be important. 

2.3 Strong optically thick wind 

In this phase, two more processes change the total angular momentum and mass of the system: a 
strong optically thick wind, and factional deposition of orbital energy into that wind. 

2.3.1 Loss of angular momentum due to an optically thick wind 

If we consider a carbon-oxygen (C-O) white dwarf accreting matter from a companion with solar 
composition, there exists a critical accretion rate above which the excess material is blown off in 
a strong wind. Hachisu et al. (1996) show that because wind velocity is about 10 times higher 




(5) 
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than orbital velocity, the wind has the same specific angular momentum as the white dwarf. The 
corresponding loss rate is 



FWIND 



q ^wind 

(1+^) Mwd 



(6) 



where q = Msg/M^d is the mass ratio and Mwind is the mass loss rate from the system. This mechanism 
will be very effective during strong optically thick wind. Note also that during the wind phase, 
the relatively small specific angular momentum content of the wind tends to stabilize mass transfer 
(Hachisu et al. 1996; Li & van den Heuvel 1997). 



2.3.2 Frictional angular momentum loss 

During an optically thick wind phase, the secondary star effectively orbits within the dense wind. Due 
to the frictional deposition of orbital energy into the wind, the separation of the components will be 
decreased. 

Livio, Govarie & Ritter (1991) estimated the change in the orbital angular momentum brought 
about by frictional angular momentum loss. 



FAML 4 U \a J Msg 



where U = (vexp/vorb) is the ratio of the expansion velocity of the envelope at the position of the sec- 
ondary to the orbital velocity of the secondary in the primary's frame of reference, a is the separation 
of the system, and Mwind is the rate of mass flow past the secondary (Warner 1995). Hachisu, Kato & 
Nomoto (1996) estimated that Vexp/vorb ~ 10- 



2.3.3 Re-accretion of material from strong wind 

We define two critical mass accretion rates onto white dwarf. First (Warner 1995): 

Mcr,i = 2.3 X 10-'^(Mwd- 0.19)3/2 [M© yr"!] , (8) 

where Mcr,i is the critical accretion rate above which hydrogen shell burning is stable. Second 
(Nomoto, Nariai & Sugimoto 1979; Hachisu et al. 1996): 

Mer,2 = 9.0 X 10-^ (Mwd - 0.5) [M© yr-1] , (9) 

where Mcr,2 is the critical accretion rate above which a strong wind occurs. In the latter case, the 
hydrogen is burning into helium at a rate close to Mcr,2, with the excess material being blown off by 
the wind at the rate: 

^wind = ^sg-^cr,2- (10) 

To calculate re-accretion of material from the wind by the secondary, we assume that, as in eq. (5), 
the mass of re-accreted material is proportional to the mass loss rate from the white dwarf (M^ind)- 

We may estimate the characteristic time scales for each of angular momentum loss mechanisms 
during a strong wind phase. We find, for example, that such a wind occurs for a binary system 
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with component masses 1.57 and 0.95 M©, and an orbital period of 1.1 d. In such a system, the 
secondary has a radius of ~2.6 R©, /i=0.4, the mass transfer rate from the secondary to white dwarf 
is 1.05x10"^ Mq yr~\ Mwind = 4 x 10~^ M© yr~\ and U=10. The characteristic time scales are 
then: Tmsw ~ 5 x 10^^ s, Tfaml ~ 3 x 10^^ s and Tfwind ~ 5 x 10^"^ s. 

3 The evolutionary code 

Models of secondary stars filling their Roche lobes were computed using a standard stellar evolution 
code based on the Henyey-type code of Paczyhski (1970), which has been adapted to low-mass stars 
(as described in detail in MS98). We use the Eggleton (1983) formula to calculate the size of the 
secondary's Roche lobe. 

For radiative transport, we use the opacity tables of Iglesias & Rogers (1996). Where the Iglesis 
& Rogers (1996) tables are incomplete, we have filled the gaps using the opacity tables of Huebner et 
al. (1977). For temperatures lower than 6000 K, we use the opacities given by Alexander & Ferguson 
(1994) and Alexander (private communication). 

3.1 Nuclear network 

Our nuclear reaction network is based on that of Kudryashov & Ergma (1980), who included the 
reactions of the CNO tri-cycle in their calculations of hydrogen and helium burning in the envelope 
of an accreting neutron star. We have included the reactions of the proton-proton (PP) chain. Hence 
we are able to follow the evolution of the elements: ^H, ^He, ^He, ^Be, ^^C, ^^C, i^N, ^^n, i^n, i^q, 
1^0, 1^0, ^^O and ^^F We assume that the abundances of ^^O and ^*^Ne stay constant throughout the 
evolution. We use the reaction rates of: Fowler, Caughlan & Zimmerman (1967, 1975), Harris at 
al. (1983), Caughlan et al. (1985), Caughlan & Fowler (1988), Bahcall & Ulrich (1988), Bahcall & 
Pinsonneault (1992), Bahcall, Pinsonneault & Wesserburg (1995) and Pols et al. (1995). 

3.2 The role of unstable helium burning 

In recent years several papers concerning unstable helium burning in accreting white dwarfs have 
been published. Cassini, Iben & Tomambe (1998) investigated accretion onto cool C-O white dwarfs 
of masses 0.5 and 0.8 M©. They considered broad range of mass accretion rates from 10~^ to 
10"^ M© yr^i. They found that even in the case when hydrogen bums at the same rate at which 
it is accreting, helium flash leads to envelope expansion. The expanded envelope, due to interaction 
with the companion star is lost from the system. The net result is that the mass of the accreting white 
dwarf decreases. 

Kato & Hachisu (1999) proposed that if an optically thick wind occurs during helium shell flashes, 
due to its high velocity (~1000 km s^^) the wind leaves the system without interaction with the com- 
panion star. They found that at a relatively high accretion rate, the mass accumulation efficiency is 
large enough for the white dwarf to grow above the Chandrasekhar mass limit. However, it is impor- 
tant to note that they considered an initially relatively massive (1.3 M©) white dwarf, and neglected 
restriction by hydrogen shell burning of the maximum accretion rate (see formula (9)). Because of 
these approximations, the mass accumulation ratio, 1 — /i, never equals unity. 
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Table 1 Initial parameters of computed sequences 



model 


Msg 




Pi(RLOF) 


X 


Z 




[M©] 


[M©] 


[d] 






A 


1.4 


0.70 


1.34 


0.68 


0.02 


B 


1.7 


0.85 


1.61 


0.68 


0.02 


C 


1.4 


0.70 


1.26 


0.50 


0.02 


D 


1.7 


0.85 


1.24 


0.50 


0.02 



Langer et al. (2000) considered the evolution of main sequence star + white dwarf binary systems 
towards type la supernovae. They assumed that above some critical accretion rate (10~^ M© yr~^ — 
for Population I) helium burning proceeds such that violent nova flashes and consequent mass ejection 
from the white dwarf are avoided. 

In our calculations, mass loss due to unstable helium burning has been neglected. So that the 
masses of white dwarfs we obtain are overestimated. 

4 Results of calculations 
4.1 The grid of models 

To understand the evolution of close binary system consisting of a C-O or NeMgO white dwarfs 
and near-main-sequence or helium-rich stars, we computed various evolutionary sequences for three 
chemical compositions (Population II — X=0.756, Z=0.001; Population I — X=0.68, Z=0.02; and 
helium rich — X=0.5, Z=0.02), for four initial mass ratios = Msg.i/Mwd.i = 1-5, 2.0, 2.5, 3.0 and 
five white dwarf masses: 0.7, 0.85, 1.0, 1.15, 1.3 Mq. For each system three different evolutionary 
stages for secondary star are followed: i) the secondary fills its Roche lobe as a M-S star with central 
hydrogen content ~ 0.4; ii) the secondary fills its Roche lobe as a terminal M-S star with small 
helium core MneCore ~ 0.01 M©; and iii) the secondary fills Roche lobe upon reaching helium core 
more MneCore ~ 0. 1 M© (Sarna, Ergma & Gerskevits, in preparation). We found that only a few evo- 
lutionary sequences are able to produce system parameters similar to U Sco. Since the evolutionary 
stage of the secondary (hydrogen- or helium-rich) is not observationally well-determined, we selected 
from our grid of models four different sequences with hydrogen and helium rich secondaries. Initial 
parameters of the sequences are presented in Table 1 . 

The sequences C and D represent the helium rich channel proposed by Hachisu et al. (1999). In 
Table 2 results for computed sequences are shown. The masses of the subgiant and the white dwarf are 
given at the moment when orbital period of the system is equal to 1.23 d. The effective temperature 
and luminosity are for the subgiant star. 

An extended grid of multi-cycle nova evolution models by Prialnik & Kovetz (1995) allows to 
find entire range of observed nova characteristics including recurrent and symbiotic novae. In Fig. 
1, three characteristic lines are drawn, according to the PK95 calculations and their classification 
scheme. Fast, slow and symbiotic novae are located below line 1 (we used the classification scheme 
of PK95). Between lines 1 and 2, recurrent novae with recurrence periods longer than 10 yr but 
shorter than 100 yr are located. Above line 2 but below line 3, the mass accretion rate onto the white 
dwarf is higher than i and nova outbursts do not occur — hydrogen burning is stable. Above line 
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Table 2 Results for computed sequences 



model 


Msg 






logL/L© 


Msg 






[M©] 


[Mq] 


[K] 




[Moyr 


A 


0.545 


0.983 


3.669 


0.153 


2.02x10" 


-8 


B 


0.936 


1.304 


3.693 


0.416 


3.58x10 


-8 


C 


0.581 


1.263 


3.719 


0.370 


4.07x10 


-8 


D 


1.696 


0.852 


3.957 


1.707 


1.14x10" 
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Table 3 Time-scales for evolutionary phases 



model 


Arni 


Afsl Afw Afs2 




Ate 






[log(A?/yr)] 






A 


4.00 


6.57 - 




8.17 


B 


5.97 


6.46 - 


7.51 


6.86 


C 


4.76 


5.41 5.39 6.17 


7.17 


6.72 


D 


4.89 


5.02 5.60 5.58 







3, the stable hydrogen burning is accompanied by a strong optically thick wind. According to our 

calculations, all models go first through an short nova phase (nl). After that, all systems enter into 
the first stable hydrogen shell burning stage (si). Sequences A and B do not exhibit a wind phase. 
After the stable burning stage, system A evolves into second nova phase (n2), while system B evolves 
into a recurrent nova phase. Sequences C and D have a wind phase (w) and a second stable hydrogen 
shell burning stage (s2). After the second stable hydrogen burning stage, sequence C evolves through 
recurrent novae phase (m). Sequence D avoids a second nova phase. For this system, the white dwarf 
mass exceeds Chandrasekhar limit during stable hydrogen burning, and a supernova explosion may 
occur. In Table 3 the duration of each phase is shown. The duration of the stable hydrogen burning 
and wind phases varies among different sequences from several hundred thousand to several million 
years. 

We infer that sequences B (hydrogen-rich) and C (helium-rich) may evolve into systems like U 
Sco. After stable hydrogen burning both sequences enter into a recurrent nova phase (Fig.l). 

For the same two evolutionary sequences, in Fig. 2 we present evolution of the orbital period 
Porb versus mass of the subgiant Msg (Fig. 2a) and mass of the white dwarf Mwd (Fig. 2b). From our 
calculations we find that both systems (B and C) evolve through the orbital period Porb = 1 -23 d twice. 
During the first crossing of Porb = 1-23 d line, the luminosity of the secondary is too high (logL/L© ~ 
0.9 and 1.35 for sequences B and C, respectively) and it does not fit the observed absolute magnitude 
Mv=+3.8 of U Sco. During the second crossing, we find from the grid of nova models (PK95 and 
KP97) outburst amplitudes A ~ 7.6 mag and recurrence periods Prec of 23 and 54 yrs, respectively, 
for sequences B and C. The mean recurrence interval implied by known outbursts is Prec = 13 yrs 
(excluding outburst in 1863). Therefore, we conclude that sequence B provides the best fit to the 
observed parameters of U Sco. 
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Figure 1: The evolution of the mass transfer rate versus white dwarf mass: (a) for a hydrogen-rich 
donor: sequence A (solid line) and sequence B (dashed line) from Table 1 ; (b) for helium-rich donor: 
sequence C (solid line) and sequence D (dashed line), also from Table 1. The lines marked 1, 2 
and 3 show lower boundaries for recurrent novae, stable hydrogen burning and strong wind regions, 
respectively. The region bounded by lines 1 and 2 shows the recurrent nova phase. The arrows show 
the direction of evolution. r> 




Figure 2: The evolution of the orbital period of the system as a function (a) of the subgiant mass Msg: 
solid line — sequence C, dashed line — sequence B, and (b) of the white dwarf mass M^d: solid line 
— sequence C, dashed line — sequence B. Horizontal thin dashed lines mark the orbital period of 
U Sco. Open circles mark positions of subgiant and white dwarf at this period for sequence B, open 
triangles for sequence C. The arrows show the direction of evolution. 
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Table 4 Chemical compositions (by mass) for sequences B and C at Porb = 1-23 d. 



model 






14n 


15n 


160 






N/C 




[xlO-3] 


[xlO-'^] 


[xlO-3] 


[xlO-"^] 


[xlO-3] 


[xlO-5] 






B subgiant 


1.41 


3.05 


3.16 


5.64 


9.77 


0.58 


4.6 


1.8 


B ejecta 


2.81 


11.51 


42.06 


251.60 


1.77 


200.00 


2.4 


10.6 


C subgiant 


0.31 


1.07 


4.69 


2.91 


9.77 


1.10 


2.9 


11.3 


C ejecta 


1.93 


7.10 


35.89 


94.82 


3.82 


51.75 


2.7 


13.6 



4.2 Chemical composition of the subgiant and ejected matter 

From the grid of models calculated by PK95 and KP97, we can interpolate the chemical composition 
of each model during its nova phase. In Table 4, the isotopic compositions in the envelope of the 
subgiant and in the ejected matter are shown (at Porb = 1-23 d) for the two sequences B and C. 

If we compare the helium-rich model C with the hydrogen-rich model B, we see that at Porb = 1- 23 
d the N/C ratio in the envelope of the subgiant is higher for the helium-rich model, but the isotopic 
ratio ^^C/^^C is higher for the hydrogen-rich model. In the ejected matter, both ratios are similar for 
models B and C. 

Our calculations show that as these models evolute, the He/H ratio of the matter lost from the 
system changes from 0.56 to 1.26. For sequences B and C, at Porb= 1-23 d the He/H ratio is 0.75 and 
0.65, respectively, i.e., within the limits (0.4 - 2) obtained from observation (Anupama & Dewangan 
2000, Williams et al. 1981). 

Figs. 3 and 4 show the evolution of the abundances by mass of helium, hydrogen, carbon, nitrogen 
and oxygen of the subgiant envelope and ejected matter. The vertical thin dashed lines show the place 
where the orbital period is equal to 1.23 d. During stable hydrogen shell burning phases, no loss of 
matter from the system occurs. These phases are marked in Fig. 4 by setting logXi = 0. 

5 Observational tests 

Chemical composition and isotopic analysis may give more information about the evolutionary stage 
of U Sco. Unfortunately, the subgiant component in U Sco is too faint for infrared spectroscopic 
observations of CO bands in order to determine the ^^C/^^C ratio. However, we think that blue 
and red domain spectra of U Sco could show some absorption structure in the region of 4216A and 
7920A in the CN sequence, similar to that observed in DQ Her (Chanan, Nelson & Margon 1978, 
Schneider & Greenstein 1979, Willimas 1983, Bauschlicher, Langhoff & Taylor 1988). Analysis of 
the blue region of the spectra is more complicated because the structure of the CN violet system can 
be affected by absorption features from the disc. However, we suggest that since the matter in the 
accretion disc reflects the chemical composition of the subgiant star, analysis of the disc will provide 
the information we need if we use observations made during the quiescent phase. We suggest that 
the red CN band is more useful for observations. We propose to observe spectral region near 7920- 
7940A to identify two ^^CN lines at 7921. 13A and 7935.67A which are important for both ^^CPc 
isotopic ratio and N/C ratio determination. In the case of low ^^C/^^C ratio (< 10), these lines are 
clearly recognized (Fujita 1985), whereas in stars with ^^C/^^C > 20 both lines are undetectable. 
Chemical analysis of the expanding envelope (Anupama & Dewangan 2000) also would give useful 



11 



n 1 1 1 1 1 1 r 

a) 



"T 1 1 1 1 1 1 1 r 







H 
He 





C 



X 
o 



-4 

-6 




n 


2 


rn 


si 




nl 













H — \ — \ — M — \ — \ — I — — \ — I — I — \ — — I — I — \ — \- 



H, He 



-4 



rn 





c 



s2 w sl.nl 



-6 



_l I I I I I u 



0.5 1 1.5 



Figure 3: The evolution of the red giant surface abundances (by mass) of H, He, C, N and O as a 
function of the subgiant mass Msg: upper panel — sequence B, lower panel — sequence C. Vertical 
thin dashed lines mark the position of subgiant mass for sequences B (Msg = 0.926 M0) and C 
(Msg = 0.603 M0) at Porb = 1-23 d, respectively (see Table 2 for more details). The first and second 
novae phases are marked nl and n2, respectively ; The first and second stable hydrogen shell burning 
episodes are marked si and s2, respectively; andjt^e strong optically thick wind phase is marked w. 
Systems evolve from right to left. For more explanation see text. 
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Figure 4: The evolution of ejecta abundances (by mass) of H, He, C, N and O as a function of the 
subgiant mass Msg: upper panel — sequence B, lower panel — sequence C. Systems evolve from 
right to left. The various phases are the same as in Fig.3 
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information, allowing for comparison with theoretical models (see Table 4). 

6 Discussion 

6.1 Hydrogen-rich vs helium-rich models 

Hachisu et al. (1999) proposed a new evolutionary path to SNe la formation, in which the companion 
star is helium-rich. In their model, typical orbital parameters of SNe la progenitors are: M^d = 
1.37 M0, Msg ~ 1.3 M©, and Msg ~ 2 x 10~^ Mq yr~^ Based on light-curve analysis, Hachisu 
et al. (2000) also constructed a detailed theoretical model for U Sco. They found that the best fit 
parameters are: M^d ~ 1-37 Mq, and Msg ~ 1.5 Mq (a range from 0.8 to 2.0 Mq being acceptable). 

However, the helium-enriched model poses a serious problem. Truran et al. (1988) discussed the 
composition dependence of thermonuclear runaway models for recurrent novae of U Sco-type. They 
showed that for M^d = 1-38 Mq, Msg = 1 .5 x 10~^ M© yr^ and L = 0. 1 L© optically bright outbursts 
are obtained only for matter with He/H<l by mass. Their results are consistent with our sequences 
B and C where He/H are equal 0.47 and 0.98, respectively. The large helium content observed in U 
Sco must then be attributed to helium enrichment from the underlying white dwarf (Prialnik & Livio 
1995). 

According to Kato (1996), a supersoft component in UV spectrum is predicted to be observable 
about 10 days after the outburst. For the hydrogen-rich model (He/H=0.1) the supersoft X-ray compo- 
nent is expected to rise until ~ 50 days after the outburst to a maximum luminosity of ~ 3 x 10^^ erg 
s^^ (J/kpc)^. For the helium-rich model (He/H = 2), the maximum luminosity is reached about 20 
days after the optical outburst. The evolution of the X-ray luminosity could therefore give important 
evidence about the chemical composition of the accreted matter. Unfortunately, only a synoptic soft 
X-ray detection of U Sco exists, some 19-20 days after outburst (Kahabka et al. 1999). 

6.2 The recurrent novae phase 

In MS98, a grid of evolutionary sequences for qi = 1.25 and white dwarf initial masses from 0.8 to 
1.2 Mq has been calculated. According to those calculations it is possible to produce a system like 
U Sco (for Population I chemical composition), but only if the initial white dwarf mass is already 
very high (> 1.2 Mq). For the model with M^da = 1-2 Mq, the system spends less than 1% of its 
semidetached evolution (5 x 10^ yr) in the recurrent novae phase. For the rest of the time, the system 
evolves through fast novae and moderately fast novae phases. Our recent calculations for qi — 1.25 
show that systems with Mwd,i = 1.3 Mq evolve towards the Chandrasekhar limit (Supemovae la). The 
system spends about 30% of its semidetached evolution (4.2 x 10^ yr) in the recurrent nova phase. 
In the calculations described here, which begin with higher initial mass ratios and lower-mass white 
dwarfs, sequences B (hydrogen-rich) and C (helium-rich) both evolve through the recurrent novae 
phase, and spend in this phase 75% and 65%of their semidetached evolution, respectively (see Table 
3). In their final phases of semidetached evolution, massive white dwarfs near the Chandrasekhar 
limit are formed, but they never exceed this limit. 
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7 Conclusions 



We calculated several evolutionary sequences to reproduce the orbital and physical parameters of the 
recurrent novae of U Sco-type. The results of these calculations can be summarized as follows: 

(1) We find a new evolutionary channel for the formation of the U Sco-type systems. We show 
that they may form from binaries with a low-mass C-0 white dwarf as accretor if Mwd,i < 0.85 M0 
and 1.5 < < 2.5. Such systems evolve through several observable stages: SSS with stable hydrogen 
burning, strong wind phase, and long-term recurrent nova phase (longer than 10^ years. Table 3). In 
the final phase of evolution, massive white dwarfs near the Chandrasekhar mass limit are formed, but 
they never exceed this limit. 

(2) We find that our evolutionary sequence B is able to produce a binary system with parame- 
ters similar to U Sco. This model has initial parameters: Mgg i = 1.7 M©, M^d,! = 0.85 M© and 
Pi(RLOF) = 1.61 d. Based on this evolutionary model, we found that the best fit parameters for U 
Sco are: Msg = 0.94 M©, M^d = 1-31 M©, logLsg/L© = 0.42, and Msg = 3.58 x 10"^ M© yr^^ for 
Porb = 1-23056 d. 
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